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ABSTRACT 

Context. SAO 244567, the exciting star of the Stingray nebula, is rapidly evolving. Previous analyses suggested that it has heated up 
from an effective temperature of about 21 kK in 1971 to over 50kK in the 1990s. Canonical post-asymptotic giant branch evolution 
suggests a relatively high mass while previous analyses indicate a low-mass star. 

Aims. A comprehensive model-atmosphere analysis of UV and optical spectra taken during 1988-2006 should reveal the detailed 
temporal evolution of its atmospheric parameters and provide explanations for the unusually fast evolution. 

Methods. Fitting line profiles from static and expanding non-LTE model atmospheres to the observed spectra allowed us to study the 
temporal change of effective temperature, surface gravity, mass-loss rate, and terminal wind velocity. In addition, we determined the 
chemical composition of the atmosphere. 

Results. We find that the central star has steadily increased its effective temperature from 38 kK in 1988 to a peak value of 60 kK 
in 2002. During the same time, the star was contracting, as concluded from an increase in surface gravity from log g = 4.8 to 6.0 and a 
drop in luminosity. Simultaneously, the mass-loss rate declined from log (M/M e yr _1 ) = -9.0 to - 1 1 .6 and the terminal wind velocity 
increased from v„ = 1800 kms -1 to 2800 kms" 1 . Since around 2002, the star stopped heating and has cooled down again to 55 kK 
by 2006. It has a largely solar surface composition with the exception of slightly subsolar carbon, phosphorus, and sulfur. The results 
are discussed by considering different evolutionary scenarios. 

Conclusions. The position of SAO 244567 in the log T e g-\ogg plane places the star in the region of sdO stars. By comparison with 
stellar-evolution calculations, we confirm that SAO 244567 must be a low-mass star (M < 0.55 M 0 ). However, the slow evolution 
of the respective stellar evolutionary models is in strong contrast to the observed fast evolution and the young planetary nebula with 
a kinematical age of only about 1000 years. We speculate that the star could be a late He-shell flash object. Alternatively, it could 
be the outcome of close-binary evolution. Then SAO 244567 would be a low-mass (0.354 M 0 ) helium pre-white dwarf after the 
common-envelope phase, during which the planetary nebula was ejected. 

Key words, stars: abundances - stars: evolution - stars: AGB and post-AGB - stars: individual: SAO 244567 - 
stars: fundamental parameters - planetary nebulae: individual: Stingray nebula (Henize 3-1357) 


1. Introduction 

The long stellar evolutionary time scales mean it is in gen- 
eral impossible for an astronomer to “watch” a star evolving 
in real time. Intermediate-mass stars ( M Z ams = 0.8-8 M 0 ) ex- 
perience their most rapid evolution close to the end of their 
nuclear-burning phase. Observing stars during this period pro- 
vides the unique opportunity to investigate on the stellar asymp- 
totic giant branch (AGB) and post-AGB evolution, including 
the AGB mass-loss phase, and the ejection, shaping, and exci- 
tation of planetary nebulae (PNe) - phases that are still not fully 
understood. 

* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26666. 

** Based on observations made with the NASA-CNES-CSA Far 
Ultraviolet Spectroscopic Explorer. 

*** Tables 1 and 3, and Figs. 4-9 are available in electronic form at 
http : //www . aanda . org 


SAO 244567, the exciting star of the Stingray nebula (Henize 
3-1357, Henize 1976), is an unusually fast evolving star. It 
was first classified to be a hot post-AGB star (Parthasarathy & 
Pottasch 1989) based on the discovery of a circumstellar dust 
shell with far-IR (IRAS) colors and flux distribution similar to 
that of PNe. Based on a spectral classification of the optical spec- 
trum obtained in 1971, Parthasarathy etal. (1995) concluded that 
the star was a B1 or B2 supergiant. From the UBV colors and 
the 1971 spectrum, they estimated T e g = 21 kK. However, they 
found that the optical spectra from 1990 and 1992 as well as the 
IUE 1 spectra (1992-1996) display many nebular emission lines, 
indicating that SAO 244567 has turned into a central star of a PN 
(CSPN) within a time span of only 20 years. Furthermore, when 
comparing the IUE spectra from 1988 and 1995 Parthasarathy 
et al. (1995) discovered it is the only known CSPN that faded 
by a factor of 2.83 in its flux level within seven years. In 
addition, it was possible for them to observe how the stellar 
wind gradually decreased. From the IUE spectrum obtained 

1 International Ultraviolet Explorer. 
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in 1988, Parthasarathy et al. (1995) measured a terminal wind 
velocity of Voo = 3500kms _1 from the Civ resonance dou- 
blet. The IUE spectrum in 1994 showed that the stellar wind 
vanished. Based on the IUE spectra, they estimated that T e g 
must be around 55 kK. Assuming a distance of 5.6kpc (Kozok 
1985b) and an expansion velocity of 8kms _1 , Parthasarathy 
et al. (1993) found that the post-AGB time of SAO 244567 is 
about 2700 years. Furthermore, they estimated the luminosity 
and core mass of the CS to be 3000 L Q and 0.55 M e , respectively. 

The first optically resolved images of the Stingray nebula 
were presented by Bobrowsky (1994) and Bobrowsky et al. 
(1998) using the Wide Field and Planetary Camera 1 (WFPC1) 
and WFPC2, respectively. Bobrowsky (1994) found that in H/3, 
the Stingray nebula appears to have an equatorial ring of en- 
hanced density tilted approximately 56° from the line of sight. 
In addition, he found bubbles of gas above and below the ring 
with areas of decreased brightness near the poles where a fast 
stellar wind has broken through the red giant envelope. From the 
H p flux, he derived an ionized mass of 0.2 A/ 0 , L = 5000 Lq 
and a stellar core mass of 0.59 M 0 . Thanks to the superior spa- 
tial resolution of the WFPC2 Bobrowsky et al. (1998) found ev- 
idence of collimated outflows, which are focused by the neb- 
ula bubbles and function like nozzles, with gas leaving through 
the polar holes. They also report a possible detection of a late 
type companion star at a distance of 2200 AU from the central 
star. 

Umana et al. (2008) present the first detailed radio study of 
the Stingray nebula by using the Australian Telescope Compact 
Array (ATCA). They find that the Stingray nebula is still em- 
bedded in the dusty remnant of the AGB phase. Depending on 
their models, they derived an ionized mass of 0.057-0.07 M Q 
and a total dust mass of 2 x 10 “ 4 M 0 in the case of silicates 
and 7.5 x 10 -5 M 0 in case of graphite. 

Arkhipova et al. (2013) performed an analysis of the neb- 
ula spectra taken in 1990, 1992, and 2011. They find significant 
changes in the relative line intensities. The low-excitation [Ol], 
[Oil], and [N II] lines became stronger relative to H/3 by a fac- 
tor of two, while the [O III] lines weakened by a factor of ~2. 
Using a formula of Kaler (1978), they estimated that T e s de- 
creased from 1990 (T eff = 57 kK) to 201 1 (T eS = 40 kK). 

Parthasarathy et al. (1995) first discovered that the observed 
properties of SAO 244567, e.g. the rapid changes in r e ff and the 
drop of luminosity, contradict with canonical post-AGB evolu- 
tion. For such a rapid evolution, the core mass should be 0.8 M Q 
or even more (Parthasarathy et al. 1995; Bobrowsky et al. 1998). 
The evolutionary time scales for CSs with core masses of 0.6 M Q 
or less are predicted to be much longer (Blocker 1995). 

To address the evolution of the properties of SAO 244567 
quantitatively for the first time, we carried out a spectral analysis 
based on all available spectra from 1988 until 2006 taken with 
IUE, FUSE 2 , HST/STIS 3 4 , and HST/FOS 4 . The comparison of 
the results to different evolutionary models should help provide 
conclusions on the nature of SAO 244567. 

This paper is organized as follows. In Sect. 2, we describe the 
observations. The spectral analysis follows in Sect. 3. In Sect. 4, 
we summarize our results and derive the distance and the mass of 
SAO 244567. We discuss possible stellar evolutionary scenarios 
and compare SAO 244567 to other low-mass CSPNe. We con- 
clude in Sect. 5. 


2 Far Ultraviolet Spectroscopic Explorer. 

3 Hubble Space Telescope/Space Telescope Imaging Spectrograph. 

4 Faint Object Spectrograph. 



2. Observations 

SAO 244567 was observed with various telescopes (Table 1). 
Parthasarathy et al. (1995) reported the decrease in the flux level 
of the IUE observations by a factor of 2.83. We found that 
SAO 244567 has faded even more; e.g., we found a decrease by 
a factor of five by comparing the IUE observation from 1988 
to the FOS observation in 1997. In Fig. 1, we show all avail- 
able, reliable flux-calibrated observations of SAO 244567. The 
STIS observations are not photometric because of the narrow slit 
width used and therefore are not shown in Fig. 1 . Comparing the 
FOS observation in 1997 to the FUSE observations in 2002, we 
found that the flux must have increased slightly. Comparing the 
flux level of the FUSE spectrum from 2002 to the one in 2006, 
we found a decrease of 15%. We mention all the available spec- 
tra and their special features in Sect. 2. 1 . The determination of 
the interstellar reddening follows in Sect. 2.2. 

2. 1. Description of the spectra 

IUE spectra were taken from 1988 to 1996. All these 
observations but SWP51772 and SWP55690 were obtained in 
low-resolution mode. The main features in the spectra are neb- 
ular emission lines (from 1992 on) and the Civ and Nv reso- 
nance doublets that are blue shifted and show P-Cygni profiles 
(Sect. 3.4). 

The HST/FOS spectra were obtained in 1997 using the 
G190H, G270H, G400H, and G570H gratings (total exposure 
time ~4600s, ProgID 6039). The observations cover a wave- 
length range from ~1600 to 6800 A. 

SAO 244567 was observed with HST/STIS in 1998 
(ProgID 7652), 1999 (ProgID 7652, 7653), and 2001 (ProgID 
8929) using the 52" x 0705 aperture. These observations cover 
a total wavelength range from 1540 to 6000 A and from 9050 to 
9650 A. In the optical wavelength range (3200 to 6800 A), the 
HST/FOS and HST/STIS spectra show many nebular emission 
lines, but also some photospheric He II lines (partially blended 
by the nebula emission lines). The NUV spectra turned out to be 
more useful for the spectral analysis, because of the He II lines 
(«—>«' = 3 — > 6, ..., 14), which could be used to determine 
log g, Teff, and the H/He ratio. We identified Cm and Civ lines 
which allowed us a more precise determination of T e ff by eval- 
uation of the Clll/ClV ionization equilibrium (see Sect. 3.2). 
Apart from the photospheric lines, we found some interstellar 
Fe n. Mg n, and Mn n lines. 

FUV spectra were taken with FUSE in 2002 and 2006 (to- 
tal exposure time: 31 ksec, ProgIDs Q304 and U109) using the 
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Fig. 2. Example for the determination of E B _ V for SAO 244567. The 
FOS observation (gray) is compared to TMAW model fluxes with dif- 
ferent corrections for the interstellar reddening. The thick (red in online 
version) line represents the TMAW model corrected with £ g _y = 0.14, 
the dashed lines (blue) indicate the error limits. 


LWRS aperture. These spectra show many interstellar lines, but 
also photospheric H I, He II, C III, C IV, N III (2002 only), N IV, 
Oiii, Oiv, Si iv, Pv, Siv, Sv, Svi, Fev, and Fevi lines. 
Furthermore, the O VI doublet lines show weak P-Cygni profiles. 
We checked the FUSE observations for relative motion of inter- 
stellar and photospheric lines, but we could not find any hint of 
velocity shifts either by comparing the observation from 2002 
and 2006, or in the individual exposures of the corresponding 
years. This indicates (except in case of a very high inclina- 
tion angle) the probability of a close companion is rather un- 
likely. All these space-based observations were retrieved from 
the MAST archive. 

2.2. Interstellar reddening 

We derived the interstellar reddening for each observation. The 
flux shape of the IUE, STIS, and FOS spectra can only be repro- 
duced using the LMC reddening law of (Howarth 1983). The 
FUSE observations could be fitted with the reddening law of 
Fitzpatrick (1999). Figure 2 shows the determination of E B ~v 
for the FOS observation. Since the impact of the interstellar red- 
dening is negligible in the infrared, the TMAW model flux of our 
best fit model to the FOS observation (T e g = 50 kK, loggr = 5.5, 
Sect. 3.2) was firstly normalized to the reddest part of the spec- 
trum. Then, we corrected the model flux for different values 
of Eb-v to fit the model flux to the observation. 

We confirm the variations in E B ~ v found by Arkhipova et al. 
(2013), who measured the reddening from the H/3 line intensity, 
using the simple relation, E B ~v = 0.68x c(H fi), with the extinc- 
tion coefficient c(H/3). In Table 2, we compare our values with 
those of other authors. 


3. Spectral analysis 

The procedure of our spectral analysis is the following. First, we 
derived T ef f, loggi, and the abundances (or at least upper lim- 
its) of the elements H, He, C, N, O, Ne, Si. P, S, Fe, and Ni 
with the Tubingen NLTE model-atmosphere package (TMAP'’, 
Werner et al. 2003; Rauch & Deetjen 2003) from the FOS, STIS, 
and FUSE spectra. The mass-loss rate of SAO 244567 during 
those years, in which these observations were taken, is negli- 
gible and therefore it is justified to use a code for hydrostatic 


5 http : //astro . uni - tuebingen . de/~TMAP 


Table 2. Interstellar reddening of SAO 244567. 


Year 

Telescope 

E B -v 

Method 

1980 

61 cm ESO 

0.18 

UBV photometry'' 

1988 

IUE 

0.18 

LMC reddening law" 

1988 


0.20 

Seaton reddening law c 

1990 

1.5 mESO 

0.18 

H /? line intensity* 



0.18 

H j8 line intensity'' 



0.14 

B aimer decrement 3 

1992 

1.5 mESO 

0.13 

H /? line intensity* 

1993 

IUE 

0.17 

LMC reddening law" 

1994 

IUE 

0.19 

LMC reddening law" 

1995 

IUE 

0.20 

LMC reddening law" 

1996 

IUE 

0.20 

LMC reddening law" 

1997 

FOS 

0.14 

LMC reddening law" 

1999 

STIS 

0.18 

LMC reddening law" 

2002 

FUSE 

0.11 

Fitzpatrick reddening law" 

2006 

FUSE 

0.11 

Fitzpatrick reddening law" 

2011 

1.9 m SAAO 

0.24 

H jS line intensity* 


Notes. Typical errors for E B -y are ±0.02. (o) This work. w Arkhipova 
et al. (2013). (c> Parthasarathy et al. (1993). rd) Kozok (1985a). 


atmospheres. We compared the PoWR 9 (Potsdam Wolf-Rayet 
model-atmosphere code, Grafener et al. 2002; Hamann & 
Grafener 2003, 2004) and TMAP models for the FUSE obser- 
vations and found that apart from the O VI P Cygni profiles, the 
spectral lines are reproduced equally well. The TMAP analysis is 
described in Sect. 3.1. To study the mass-loss rates and terminal 
wind velocities from the P Cygni profiles found in the IUE and 
FUSE spectra and to determine T e ff in the years of the IUE ob- 
servations, we used PoWR. The description of this analysis can 
be found in Sect. 3.4. 

3.1. TMAP model atmospheres 

We used TMAP to compute non-LTE, plane-parallel, fully 
metal-line blanketed model atmospheres in radiative and hydro- 
static equilibrium. The final models included opacities of the el- 
ements H, He, C, N, O, Ne, Si, P, S, Fe, and Ni. The model 
atoms for this analysis were taken from the Tubingen model- 
atom database TMAD and calculated (Fe, Ni) via the Tubingen 
iron-group opacity interface TIRO 8 (Muller- Ringat 2013). It 
has been developed recently in the framework of the Virtual 
Observatory (VO 9 ) and is provided as a registered service by 
the German Astrophysical Virtual Observatory (GAVO 10 ). The 
statistics of the model atoms used in this analysis are summa- 
rized in Table 3. All spectral energy distributions that were calcu- 
lated for this analysis are available via the registered GAVO ser- 
vice TheoSSA * 11 . 

In a standard procedure (e.g. Ziegler et al. 2012), we mod- 
eled both, the photospheric and the interstellar line-absorption 
spectrum, to correctly identify the stellar lines. We employed the 
OWENS program to model the ISM line absorption with sev- 
eral ISM clouds with different parameters (radial and turbulent 
velocities, temperatures, and column densities of the individual 
ions). In the following, we describe our photospheric analysis. 


6 http : //www. astro .physik. uni -potsdam. de/~wrh/PoWR/ 

7 http : //astro . uni - tuebingen . de/~TMAD 

8 http : //astro . uni - tuebingen . de/~TIR0 

9 http://www.ivoa.net 

10 http://www.g-vo.org 

11 http://dc.g-vo.org/theossa 
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Fig. 3. Comparison of our best-fit TMAP models (black, red in online version) to the FOS and STIS observation (gray). Hell, ClII-IV, and 
interstellar Fe II lines are marked. The vertical bar indicates 10% of the continuum flux. 


3.2. Effective temperature, surface gravity, H/He ratio, 

C abundance 

Based on the IUE spectrum in 1995, Parthasarathy et al. (1995) 
estimated T e ff ~ 55 kK. Starting with this value, we computed a 
grid of H+He composed model atmospheres (7" e ft = 40-70 kK, 
10 kK steps), log g = 4.0-7. 0, 0.5 steps), H/He abundance ratios 
(by mass) of 9/1, 3/1 (i.e. solar), and 2/3 for a preliminary T e s 
and log g determination. The synthetic spectra of these models 
were then compared to the observed He line profiles in the STIS 
and FOS spectra. 

Hen A 2723.2 A is much too strong in all H/He = 2/3 mod- 
els so we can rule out He enrichment. The H/He = 9/1 models 
show that SAO 244567 is not He-poor because it is not pos- 
sible to match the strengths and wings of the Hell lines. All 
line profiles are reproduced well with a solar H/He ratio at 
log g = 6.0 and T e u = 50, 60, 70 kK for the STIS spectrum. For 
the FOS observations the best fit was found for \ogg = 5.5 and 
T e ff = 40, 50, 60 kK. We adopted this abundance ratio for our 
analysis. 

We identified several Cm and Civ lines in the FUSE, 
STIS and FOS spectra. Consequently, we included C in our 
models to evaluate the C Ill/C IV ionization equilibrium and to 
determine the C abundance. For the determination of the C 
abundance, we used Cm AA 1 165.6-1 165.9, 1 174.9-1 176.4 A 
and Civ AA 1107.6-1108.0, 1 168.8-1 169.0 A in the FUSE 
spectrum and Chi 42297. 6A and Civ AA 2336.7-2337 .1, 
2405.1-2405.9, 2524.41-2530.7 A in the STIS and FOS spec- 
tra. The C iv AA 1548.2, 1550.8 A resonance doublet is contam- 
inated by the respective interstellar lines, hence not suited to 
determine the C abundance. The best agreement with the ob- 
servation was found at C = 5.0 x 10 , which is 0.2x solar 
(solar values according to Asplund et al. 2009) in all spectra. 
Simultaneously, we refined our model grid (T e ft step 5kK). In 
Fig. 3, we show the best-fit models for the FOS observation in 
1997 and the STIS observation in 1999. We discovered that the 
ionization equilibrium and also the surface gravity have changed 
with time. While we determined T e s = 50 kK and log g = 5.5 for 


Table 4. Temporal evolution of the atmospheric parameters as derived 
from our spectral analysis. 


Year 

T t ff 

[kK] 

lo g g 
[cm/s 2 ] 

log M 
[M 0 /yr] 

Voo 

[km s -1 ] 

Code 

1988 

38 

4.8 

-9.0 

1800 

PoWR 

1992 

43 

5.0 

-9.0 

1800 

PoWR 

1993 

44 

5.0 

-9.1 

2100 

PoWR 

1994 

48 

5.2 

-10.0 

2100 

PoWR 

1995 

50 

5.2 

-10.0 

2400 

PoWR 

1996 

50 

5.2 

-10.0 

2800 

PoWR 

1997 

50 

5.5 



TMAP 

1999 

55 

6.0 



TMAP 

2002 

60 

6.0 

-11.3 

2800 

TMAP/PoWR 

2006 

55 

6.0 

-11.6 

2800 

TMAP/PoWR 


Notes. Typical errors are Ar eff = ±5 kK, Alog g = ±0.5, AM = +0.2 dex, 
and Au„ = +200 km s -1 . 


the FOS observation, the STIS observation is best reproduced by 
Teg = 55 kK and log g = 6.0. 

In the FUSE spectra, we evaluated ionization equilibria 
of Cm/Civ, Nlll/Niv (FUSE 2002 only), Om/Oiv, and 
S IV/S v/S VI to derive Teff. In Fig. 4, we show the different ion- 
ization equilibria of the FUSE observation in 2002. The model 
with T e « = 60 kK matches the observations. In the 55 kK model, 
the O III and S IV lines are too strong, whereas the S VI line is too 
weak. In the 65 kK model, the N III, O III, S IV, and C III lines 
are too weak and the S VI line is too strong. Comparing the 
FUSE observations from 2002 and 2006, we found that the O III 
and S IV lines became slightly stronger. Figure 5 shows that the 
best agreement in 2006 is obtained for T e ff = 55 kK. The value 
of log g = 6.0 ± 0.5 was confirmed using the He II line wings. In 
Table 4, we list the values of Teff and log g that we found from 
this analysis. 

3.3. Element abundances 

For determining the abundances of N, O, Si, P, S, and Ni, we 
used the respective lines found in the FUSE spectra. To derive 
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upper limits, we compared our models to FOS spectra. For 
the determination of the element abundances, we adopted the 
values of T’eif and log g found for each epoch (Sect. 3.2). The 
Fie and C abundances were already achieved as a byproduct 
of the T eff and loggr determination (see above). To calculate 
the SAO 244567 model grids in a reasonable time, we only in- 
cluded H, Fie, C, N, and O and added the trace elements Ne, 
Mg, Si, P, S, Fe, and Ni in a subsequent line-formation calcula- 
tion; i.e., we kept the atmospheric structure fixed and calculated 
NLTE occupation numbers only for newly introduced species. 
This is justified by the comparison of temperature structures of 
a FfHeCNO and a complete, final HHeCNONeSiPSFeNi model, 
with T eff = 60 kK and log g = 6.0, where the atmospheric struc- 
ture is calculated with all these elements. Figure 6 demonstrates 
that the temperature structures are almost identical in the line- 
formation region (-4 ^ logm SI). All abundances found from 
line-formation calculations are verified with our final models. 

Nitrogen and oxygen. The N abundance of SAO 244567 
is solar (N = 6.9 x 10 -4 ). It was measured us- 

ing NmTT 1183.0-1 184.6 A, Niv AA 1123.5-1 133.1 A, and 
N iv A 1 183.0 A lines in the FUSE spectrum. The O abundance 
was derived from the Om AA 1149.6-1 153.8 A, and OlV AA 
1045.4-1050.5, 1167.8-1168.0A lines identified in the FUSE 
spectrum (Fig. 4). We found a solar value (O = 5.6 x 10 -3 ). 

Neon, sodium, magnesium, silicon, phosphor, and sulfur. 
The only Ne line that was prominent in our models is 
Ne hi A 2678.7 A. It would emerge from the STIS and FOS spec- 
tra above an upper abundance limit of five times the solar value 
(Ne £ 6.2 x HU 3 ). Using Nam TT2012.5,2031.8A, we de- 
rive an upper limit for sodium of ten times the solar value 
(Na S 2.9 x 10 -4 ). For magnesium, we derive an upper limit 
(solar. Mg < 7.1 x 10“ 4 ) using Mgm U12134.7,2178.4A. For 
the determining the Si abundance we used the Sim A 1 113.2 A 
and Siiv A 1128.3 A lines. All other Si lines are contami- 
nated by interstellar absorption lines. As illustrated in Fig. 7, 
the best fit is obtained at a solar Si abundance (Si = 
6.7 x 10 -4 ). For lower/higher Si abundances the modeled Si 
lines are too weak/strong. Pv AA 1 118.0, 1 128.0 A were used 
to determine the P abundance (Fig. 8) which is 0. 1 x so- 
lar (P = 5.8 x 10“ 7 ). We identified S iv AA 1073.0, 1073.5 A, 
S v A 1039.9 A, and S vi A 1 117.8 A and determined the S abun- 
dance to be 0.5x solar (S = 1.6 x 10 -4 ). 

Iron and nickel. Some Fe lines could be identified in the FUSE 
spectra (Fev AA 999.6-1001.1 A, Fevi AA 1165.7, 1167.7 A). 
The Fe abundance was found to be solar (Fe ^ 1.3 x 10 -3 ). 
The Fe v/Fe vi ionization equilibrium also confirms T e ff = 60 kK 
(Fig. 9). The quality of the available FUSE spectra is not suffi- 
cient to unambiguously identify individual Ni lines. Therefore, 
we could only derive a (solar) upper limit for the Ni (Ni :£ 
7.1 x 10 -5 ) abundance (Fig. 9). We summarize the abundances 
in Table 5. Their error limits are ±0.3 dex. 

3.4. Stellar wind 

Parthasarathy et al. (1993) found that in the IUE spectra of 1988 
and 1992 the Civ and N V resonance doublets are blue-shifted 
and they show P-Cygni profiles. Fooking more carefully, we 
found that all the spectra show P-Cygni profiles in these lines. 
While in 1988 the P-Cygni profile of N V is quite weak and 
C IV strong, this changes from 1992 on, which indicates that the 


Table 5. Photospheric abundances of SAO 244567. 


Element 

log mass 
fraction 

log number 
fraction 

log 6 

[Al 

H 

-0.13 

-0.04 

12.01 

0.00 

He 

-0.60 

-1.11 

10.94 

0.00 

C 

-3.30 

-4.28 

7.77 

-0.68 

N 

-3.16 

-4.21 

7.84 

0.00 

O 

-2.25 

-3.35 

8.70 

0.00 

Ne 

<-2.21 

<-3.41 

<8.64 

<0.70 

Na 

<-3.53 

<-4.80 

<7.26 

<1.00 

Mg 

<-3.15 

<-4.44 

<7.61 

<0.00 

Si 

-3.18 

-4.53 

7.52 

0.00 

P 

-6.24 

-7.63 

4.42 

-1.00 

S 

-3.81 

-5.22 

6.83 

-0.30 

Fe 

-2.89 

-4.54 

7.51 

0.00 

Ni 

<-4.15 

<-5.82 

<6.23 

<0.00 


Notes. [A] denotes log (abundance/solar abundance), log ex are normal- 
ized to log XfrxSx = 12.15, where fix is the atomic weight. 


star’s T e ff has increased. Furthermore, we found in the FUSE 
spectra that the Ovidd 1031.9, 1037. 6A resonance lines show 
P-Cygni profiles, too. 

To measure the mass-loss rate and Teff of SAO 244567, we 
used PoWR. It solves the NFTE radiative transfer in a spheri- 
cally expanding atmosphere simultaneously with the statistical- 
equilibrium equations and accounts at the same time for energy 
conservation. Iron-group line blanking is treated by means of the 
superlevel approach (Grafener et al. 2002), and a wind clumping 
in first-order approximation is taken into account (Hamann & 
Grafener 2004). We do not calculate hydrodynamically consis- 
tent models, but assume a velocity field following a /Maw with 
P = 1 . For the PoWR models, we adopted the elemental abun- 
dances that we found in our TMAP analysis. We extrapolated 
the mass and the luminosity of SAO 244567 using the evolution- 
ary tracks from Hall et al. (2013) and found M = 0.36 M 0 and 
log (T/T 0 ) = 2.5. We also calculated models with M = 0.47 M Q 
and log (L/L q ) = 2.2 (values extrapolated from the tracks by 
Driebe et al. 1998 ), but we found that the line profiles change 
only slightly with M and L. This agrees with Herald & Bianchi 
(2007) who report that wind features are not very sensitive to 
the gravity, which is in the case of the PoWR code an equivalent 
input parameter. 

Based on these assumptions, we varied the terminal wind 
velocity v mass-loss rate M, and Teff in our models. The rela- 
tive strengths of N v/C IV turned out to be very sensitive to T e ff, 
so that we can achieve an error of only AT e ff = ±3kK. The 
sensitvity to T e ff can be clearly seen by comparing the best-fit 
models for 1988 and 1992, where the only difference of these 
models is T ef j. In the T e g = 38 kK model for 1988, Civ is rel- 
atively stronger than N V, whereas the opposite holds for the 
T e ff = 43 kK model for 1992. The mass-loss rate was derived 
from the strengths of the P-Cygni profiles. As the P-Cygni pro- 
files broaden as Voo increases, Voo can be measured not only from 
the blue edge of the absorption component of the P-Cygni pro- 
files, but also from the shape of the emission peak. 

The best-fit models are shown in Figs. 10 (IUE spectra) 
and 11 (FUSE). The P-Cygni profile of Ovid 1037.6 A in the 
PoWR models is blended by interstellar Hi and is therefore not 
visible in the observation. Furthermore, we found that the ab- 
sorption components of O vi AA 1031.9, 1037. 6 A must have an 
interstellar origin. Table 4 lists the parameters that we found 
from this analysis. The resulting surface gravities are also given. 
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Fig. 10. Comparison of the IUE spectra taken from 1988 to 1996 (thin, 
gray) with the best-fit PoWR model (thick, red). All but the 1988 obser- 
vation are shifted for clarity. 



Fig. 11. Comparison the FUSE spectra taken in 2002 and 2006 (gray) 
with the best-fit PoWR model (thick, red). The blue line (thin) indicates 
the combined stellar and ISM spectrum. 

Table 6. Comparison of the element abundances (number fractions rel- 
atively to H) of SAO 244567 (CS) as derived by our analysis and of its 
PN as derived by analysis of the spectrum in 1990 (Parthasarathy et al. 
1993), 1992, and 2011 (Arkhipova et al. 2013). 


We point out that these values refer to the models with M = 
0.36 M e , the resulting values for the surface gravity from the 
M = 0.47 M q models were about 0. 1-0.2 dex higher. 

We find that the mass-loss rate decreased continuously from 
log(M/M Q yr _I ) = -9.0 in 1988 to \og(M/M 0 yr” 1 ) = -11.6 
in 2006. The previously published value (ttoo = 3500kms _1 , 
Parthasarathy et al. 1995) turned out to be an overestimate. We 
find that the terminal wind velocity has steadily increased from 
Voo = 1800 kms -1 in 1988 to Voo = 2800 kms -1 in 2006. 

4. Results and discussion 

We analyzed the high-resolution UV spectra taken with FUSE, 
STIS, FOS, and IUE by means of non-LTE line blanketed model 
atmospheres. The values of T e ff were derived with high preci- 
sion by evaluating the ionization equilibria of Cm/Civ (FOS, 
STIS and FUSE spectra), O III/O IV, S IV/S V/S VI (FUSE spec- 
tra), and N III/N IV (FUSE 2002 spectra). From the IUE spectra, 
we derived T e g using the relative strength of N v/C IV. 

Both the surface gravity and the H/He ratio (solar) were ob- 
tained by an analysis of the Hen lines found in the STIS and 
FOS spectra. For the FOS observations, the best fit was found 
with \ogg = 5.5, while for the STIS observation the model with 
log g = 6.0 fits better. This value of log g = 6.0 was confirmed 
by the wings of the He n lines in the FUSE spectra. 

SAO 244567 exhibits solar abundances of H, He, O, and N. 
The C abundance is definitely subsolar. This indicates that the 
AGB phase of the star was terminated before the third dredge- 
up (Mello et al. 2012). P and S were found to be subsolar, the 
Si and Fe abundances are solar. For Ne, Mg, and Ni, we could 
only derive upper limits (five times solar for Ne, ten times solar 


CS PN 

2011 1992 1990 


He/H 

8.50 

X 

10- 

-2 

9.60 

X 

io- 

-2 

9.30 

X 

10- 

-2 

1.03 

X 

io - 1 

C/H 

5.69 

X 

10 

-5 





7.59 

X 

10' 

-5 




N/H 

6.69 

X 

10 

-5 

5.75 

X 

io- 

-5 

6.45 

X 

10' 

-5 

6.50 

X 

10“ 5 

O/H 

4.84 

X 

10 

-4 

1.91 

X 

io- 

-4 

2.46 

X 

10' 

-4 

3.00 

X 

10“ 4 

Ne/H 

<3.59 

X 

10 

-4 

3.47 

X 

io- 

-5 

5.75 

X 

10' 

-5 

9.20 

X 

10“ 5 

S/H 

6.03 

X 

10 

-6 

2.40 

X 

io- 

-6 

6.03 

X 

10' 

-6 

2.20 

X 

10“ 6 


for Na, solar for the rest). We could not find any hint of a change 
in the chemical abundances during the years. Within the error 
limits (typically ±0.3 dex) the abundances found for the CS are 
in good agreement with the nebula abundances (Table 6) deter- 
mined by Arkhipova et al. (2013) and Parthasarathy et al. (1993). 

The temporal evolution of the atmospheric parameters of 
SAO 244567 is summarized in Table 4 and Fig. 12. We 
found that compared to the temperature found from the spec- 
trum in 1971 (Sect. 1) SAO 244567 has increased its T e ff by 
about 40 kK within only thirty years. The peak T e g of 60 kK 
was reached in 2002, and the 2006 observations suggest that it 
is now decreasing. We cannot confirm the values for T e ff found 
by Arkhipova et al. (2013). They find the peak T e ff in 1990 
(Teff = 57 kK), i.e. about 10 kK higher than our analysis from 
spectra taken at that time. We want to point out, that the formal- 
ism of Kaler (1978) used by Arkhipova et al. (2013) might not 
provide good estimates for the CS temperature. Comparing the 
T e ff of Kaler (1978) found for his sample of CSPNe, shows of- 
ten big deviations from recent literature values (e.g. NGC 6543: 
Kaler 1978 found T t ff = 45 kK while Herald & Bianchi 2011 
found Teff = 60 kK). However, evaluation of the ionization equi- 
libria found in the FUSE spectra shows J e ff declining from 2002 
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1970 1980 1990 2000 2010 

year 


Fig. 12. Temporal evolution of T t g, logg, mass-loss rate M. and the 
terminal wind velocities v m as derived in our analysis (red). The black 
stars in the bottom panel show the T e fr as derived by Arkhipova et al. 
(2013), the black square indicates the T e g estimate when SAO 244567 
was still a B-type supergiant (1971). 


on, supporting the contention by Arkhipova et al. (20 1 3) that T t g 
has recently peaked and has begun to decline. 

By examining the P-Cygni profiles of the C IV and N V res- 
onance doublets in the IUE spectra from 1988 to 1996 and the 
P-Cygni profiles of Ovi in the FUSE spectra, we derived the 
mass-loss rates and terminal wind velocities as a function of 
time. We found a steady decrease in the mass-loss rate (from 
log(M/M 0 yr _1 ) = -9.0 in 1988 to log(M/M 0 yr” 1 ) = -11.6 
in 2006) and the steady increase in the terminal wind velocity 
(from Voo = 1800kms _1 in 1988 to = 2800kms _I in 2006). 
To compare the mass-loss rate with predictions from radiative- 
driven wind theory, we used Eq. (5) in Vink & Cassisi (2002) 
with L = 2.5 L 0 , M = 0.354 M 0 , T eff = 38 kK, and Z = 2 to de- 
rive the theoretical value of the mass-loss rate in 1988. We found 
log M = -8.8 ± 0.3, which is in good agreement with the ob- 
served value. However, we note that the formula given by Vink 
& Cassisi (2002) might not give the right predictions, since it is 
meant for a slightly different parameter space in L, Teff, and M. 

The decrease in flux with a simultaneous increase of T e g 
from 1988 to 1999 can be explained by the contraction of the 
star ( L oc R 2 T e ff 4 oc g _1 T e fj 4 ). This is consistent with the increase 
in loggi as seen in the FOS (1997) and STIS (1999) observations, 
and consistent with the increase in u<x, which is proportional to 
u eS c K Q- The decrease in flux from 2002 to 2006 can be ex- 
plained by a decrease in T e ff as seen in the FUSE spectra. 



Fig. 13. Evolution of SAO 244567 (black dots) in the log T eff -log g 
plane compared to post-AGB (red, thick) by Blocker (1995), post-EHB 
(blue, thin) by Dorman et al. (1993), and post-RGB (green, dashed) evo- 
lutionary tracks by Hall et al. (2013) The tracks are labeled with stellar 
masses. The ellipse indicates the errors of T t g and log g in 1997. 

4.1. Distance 

The distance of SAO 244567 was estimated first by Kozok 
(1985b) using UBV photometry and the assumption of a typi- 
cal absolute magnitude of a Be star, leading to d = 5.64kpc. 
Arkhipova et al. (2013) argues that Kozok (1985b) overes- 
timated the absolute magnitude of SAO 244567 since it was 
not a main sequence star at that time. Using the estimates of 
Parthasarathy et al. (1993) for T e g = 37.5 kK and log g = 4.0, 
they estimated M ~ 0.55 M 0 and derived d = 1.8 kpc. We cal- 
culated the spectroscopic distance of SAO 244567 using the flux 
calibration of Heber et al. ( 1984) for d e ff = 5454 A, 

d[ pc] = 7.11 x 10 4 ■ -y/ H v -Mx io 0 - 4 '"v 0 -iog 0 , 

with my 0 = tny - 2.175c, c = \A1Eb-v , and the Eddington 
flux H v (7.15 x 10 -4 erg/cm 2 /s/Hz) at 5454A of our model at- 
mosphere (T e ff = 50 kK and loggr = 5.5). The visual brightness 
(my = 14.9) was extracted from the FOS spectrum, so we 
used all the model parameters found for this year. We assumed 
Eb-v = 0. 14 (Sect. 2.2) and a stellar mass of M = 0.354^-q 4 M & 
(Sect. 4.2). We derived d = 1.6^° * kpc which leads to a height 
above the Galactic plane of z = 0.3 ± 0.2 kpc. These values 
are in good agreement with the estimates of Arkhipova et al. 
(2013). From the angular diameter of 273, a linear radius of 
the PN of R = 0.009^5'qq 4 pc results. With the expansion ve- 
locity of 8.4km s -1 (Arkhipova et al. 2013), the kinematic age 
of the PN is only about 1013^3 years. We note that different 
masses, which we derive from various evolutionary scenarios, 
hardly affect the distance determination and that the large errors 
are mainly due to the uncertainties in log g. 

4.2. Mass and evolutionary status 

In Fig. 13, we compare the time-dependent location of 
SAO 244567 in the log T ef f-loggr plane with H-rich post-AGB 
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evolutionary tracks by Blocker (1995), post-EHB-tracks by 
Dorman et al. (1993), and post-RGB-tracks by Hall et al. (2013). 
The last represent the evolution of post-common envelope rem- 
nants of post-RGB stars that will not ignite He and become low- 
mass white dwarfs with a He core. The determination of the mass 
of SAO 244567 would only be possible, if the evolutionary status 
of this object is known. For that reason, we give only the mass 
estimates valid for a certain evolutionary scenario below. 

The position of SAO 244567 in the log Teff-log g plane 
places the star in the region of sdO stars. Since SAO 244567 does 
not match the 0.546 M G track within the error ranges, it may ei- 
ther be a very low-mass post-AGB star with a mass of about 
0.53 M G or even an AGB-manque star, which never reaches the 
AGB. 

The existence of the PN around SAO 244567 suggests that it 
is more likely a post-AGB star than an AGB-manque star be- 
cause - in general - only these stars are expected to eject a 
PN. If we consider SAO 244567 as a 0.53 M 0 post-AGB star, it 
evolves much more quickly than predicted by theory. Comparing 
T e ff = 21 kK estimated for 1 97 1 to the peak value of T t $ = 60 kK 
in 2002, SAO 244567 increased its temperature to 40 kK within 
only 30 years. According to Fig. 6 in Blocker (1995) such a 
rapid evolution (logfdt/dTeff/yrK -1 ) = -3) would be expected 
only for a 0.87 M G star with an initial mass of about 6 M Q . For 
a 0.55 M G star it would take 130000 years to heat up from 
Tea =21kKto7’ eff = 60kK. 

An alternative to the canonical post-AGB evolution is a 
late thermal pulse (FTP, e.g. Blocker 2001). The evolutionary 
speed of objects, which are considered to have undergone a 
(very) late He-shell flash (V605 Aql, e.g., Clayton et al. 2006; 
V4334 Sgr, e.g., Hajduk et al. 2005; and FG Sge, e.g., Jeffery 
& Schonberner 2006), is very high (decades). A very late ther- 
mal pulse can be ruled out because this scenario produces a 
hydrogen-free stellar surface, whereas a normal surface compo- 
sition is typical of a LTP in the phase relevant to SAO 244567. 
The FTP occurs when the star evolves with roughly constant lu- 
minosity from the AGB towards the white dwarf domain. The 
convective shell triggered by excessive helium burning is not 
able to penetrate the hydrogen-rich envelope from below be- 
cause the entropy jump across the helium/hydrogen interface 
is too large. Only when the star evolves back to its Hayashi 
limit on the AGB (T e ff < 7000 K), envelope convection sets in 
again (Blocker & Schonberner 1996, 1997; Schonberner 2008). 
Therefore, the surface abundances of SAO 244567 are in agree 
with the FTP scenario. The evolutionary calculations of Blocker 
(2001) for a 0.625 M G FTP star predict that the solar surface 
composition of the star holds up after the FTP. Only when the 
star has cooled down to T e ff =10 kK, the H abundance decreases 
while the He and C abundances increase. The FTP scenario is 
able to explain the rapid evolution (increase in T e ff and drop 
in luminosity). It would predict a decrease in T e g (the FUSE 
observation in 2006 already suggests that) and an increase in 
the brightness within the next decades (similar to the evolu- 
tion of FG Sge). The FTP is a good candidate for explaining 
the evolutionary status of SAO 244567; however, concrete evo- 
lutionary models that match the position of SAO 244567 in the 
log Teff— log g plane, are missing to prove this scenario. 

Considering SAO 244567 as an AGB-manque star, which ig- 
nited central He-burning, we derive a mass of 0.49 M G according 
to the tracks of Dorman et al. (1993). Within this scenario the 
existence of the PN and the short evolutionary time scales are 
even more difficult to explain. Even if the RGB precursor would 
have ejected a PN, it should have disappeared a long time ago 
since the horizontal branch (HB) evolution lasts ~10 8 years. It is 


worthwhile mentioning that the mass of the remaining H layer 
(<0.001 M G ) of extended horizontal branch (EHB) stars is much 
too low to produce a nebula at the end of the HB stage. We can 
also rule out an early hot-flasher scenario. Brown et al. (2001) 
found that as the Reimers mass-loss parameter g r increases, the 
peak of the main helium core flash shifts to higher temperatures, 
and the subsequent zero age horizontal branch (ZAHB) position 
becomes hotter. However their models never produced stars with 
a ZAHB position hotter than T e s =31.5 kK. 

Recently Hall et al. (2013) have proposed that it is possi- 
ble that at least some PNe are composed of matter ejected from 
a binary star system during common envelope (CE) evolution. 
For these PNe, the ionizing component is the hot and luminous 
remnant of a giant that had its envelope ejected by a compan- 
ion in the process of spiraling in to its current short-period orbit. 
A large fraction of CE phases that end with ejection of the en- 
velope are thought to be initiated by low-mass red giants, gi- 
ants with inert, degenerate helium cores. In their calculations 
Hall et al. (2013) find, that PNe are expected in post-CE sys- 
tems with core masses greater than about 0.3 M G if remnants end 
the CE phase in thermal equilibrium. Considering this scenario 
for SAO 244567, we estimate a mass of 0.354 ± 0.03 M G . The 
rapid evolution of SAO 244567, though, is also not predicted in 
their models. A higher mass-loss rate might increase the evolu- 
tionary speed, but there are no hints of any enhanced mass loss 
in the observations. However, if the remnant is in thermal non- 
equilibrium after the CE ejection, it might evolve rapidly enough 
(Philip Hall, priv. comm.). 

The current decrease in T e fi of SAO 244567 from 2002 on 
found in our work and also by Arkhipova et al. (2013) argues for 
post-RGB evolution instead of the post-EHB evolution, since the 
H-shell-burning models cool down directly, whereas the He-core 
burning post-EHB models become even hotter (Fig. 13). 

The bipolar shape of the PN is most likely due to some 
kind of binary interaction of a RGB or AGB star with either 
a stellar or substellar companion. This was already suggested 
by Bobrowsky (1994) to explain the axisymmetric structure of 
SAO 244567. Soker (1997) proposed that bipolar structures of 
PNe, as seen in SAO 244567, can be explained by close binary 
interactions that avoid a CE phase, or they entered this phase 
only in their late evolution. The gravitational interaction between 
the companion and the PN progenitor and/or their winds are non- 
negligible. The companion’s tidal force can spin up the primary 
and enhance the mass-loss rate or the companion can accrete 
mass from the primary’s wind. This accretion process is likely 
to result in a high mass-loss rate in the equatorial plane, leading 
to a bipolar PN. The possible companion detected by Bobrowsky 
et al. (1998), however, is too far away (638AU at a distance 
of 1.6 kpc) for significant interaction with SAO 244567. But it 
might be possible that SAO 244567 has another, hitherto unde- 
tected companion. 

4.3. Comparison with other low-mass CSPNe 

SAO 244567 is not the only low-mass CSPN. Napiwotzki (1999) 
found five such candidates (HDW 1 1 , K 2-2, GD 56 1 , PHF 932, 
and DeHt5) that are not in accordance with post-AGB evolu- 
tion. For the DA-type WD DeHt5 (M = 0.414 M G ), they con- 
sidered a post-RGB evolution as the kinematical age of the PN 
(~ 129 000 yrs) is in good agreement with the post-RGB age 
(~100000yrs) estimated from the calculations of Driebe et al. 
(1998). The two DAO-type WDs K2-2 (M = 0.39 M G ) and 
HDW 11 (M = 0.38 M g ) are also candidates for this evolution- 
ary scenario. 
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For PHL 932 and EGB 5, an object with similar parameters, 
the post-RGB times corresponding to the tracks of Driebe et al. 
(1998) are much longer than the ages of the PNe. The PN EGB 5 
is believed to be a remnant of an ejected CE. Geier et al. (2011) 
discovered a close low-mass companion orbiting the sdB cen- 
tral star and thus EGB 5 is considered to be a good candidate 
for a post-CE object. Napiwotzki (1999) considered PHL 932 
{M = 0.28 M q ) to be the outcome of a CE event of a star with 
a degenerate CO core as calculated by Iben & Tutukov (1985) 
or even as the result of a merger within a common envelope 
(Mendez et al. 1988b). However Frew et al. (2010) convincingly 
demonstrate that the nebula around PHL 932 is not a PN, but 
rather a Stromgren sphere in the ambient ISM. Also EGB 5 has 
recently been classified as an ionized Hn region (Frew et al. 
2013). The same holds for the nebulae around the DAO-type 
WDs GD 561 and BD-22°3467 (Ziegler et al. 2012), which Frew 
& Parker (2010) also consider to be Stromgren spheres. Other 
objects of this type are the sdO stars HD 497898, LSS630, 
LSE44, LSE153, SB 705, and KPD 0720-0003 around which 
nebulosities were detected that are not believed to be ejected by 
or physically associated with the corresponding stars (Mendez 
et al. 1988a). This can clearly not be true for SAO 244567 owing 
to its complex nebula structure and the expansion velocity. 

A multi-shell planetary nebula around the hot sdO star 
2MASS J193 10888+4324577 has recently been detected 
(Aller et al. 2013). The morphology of the nebula (bipolar and 
elliptical shell, whose major axes are oriented perpendicular 
to each other) strongly resembles the one of the Stingray 
nebula. The complex nebula structure and the fact that 2MASS 
J193 10888+4324577 was found to be in a binary system 
(Jacoby et al. 2012) suggest that it might have formed through 
binary star evolution. 


5. Conclusions 

SAO 244567 is a rapidly evolving object. Its evolutionary sta- 
tus remains unclear. The most reasonable explanations are a late 
He-shell flash or CE evolution with a remnant that is in ther- 
mal non-equilibrium after the CE ejection. However, respective 
models are lacking that match the position of SAO 244567 in the 
log Teff-log g plane. The contradiction between observations and 
theory make SAO 244567 particularly interesting. Its fast evolu- 
tion gives us the unique opportunity to study stellar evolution 
in real time and establishes constraints for stellar evolutionary 
theory. Further observations, in the next years, decades and even 
centuries, are essential for monitoring whether the rapid evolu- 
tion of SAO 244567 is still going on and to see if it is directly 
evolving to the white dwarf domain or back to the AGB. The de- 
tection of a close binary would support the scenario of a CE ejec- 
tion, whereas an increase in brightness and decrease in T e ff over 
the next decades, would indicate an evolution back to the AGB 
and hence speak for a LTP scenario. 
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Table 1 . Observations of SAO 244567. 


Instrument/ 

telescope 

Obsld 


Aperture/ 

grating 

Start time (UT) 

Exposure time (s) 

B&C, ESO 1.5 m 




1990-06 

1-600 

B&C, ESO 1.5 m 




1992-08 

1-600 

IUE 

SWP33954LL 


large 

1988-07-21 23:44:32 

1800 


LWP13715LL 


large 

1988-07-22 00:34:19 

1200 


SWP33955LL 


large 

1988-07-22 01:05:33 

1800 


LWP22874LL 


large 

1992-04-23 04:19:56 

600 


SWP44459LL 


large 

1992-04-23 03:52:00 

1200 


LWP25397LL 


large 

1993-04-23 05:06:09 

600 


SWP47530LL 


large 

1993-04-23 05:23:38 

1200 


SWP50589LL 


large 

1994-04-19 08:13:06 

1200 


LWP2797 ILL 


large 

1994-04-23 01:56:53 

600 


SWP51772HL 


large 

1994-08-11 01:52:34 

25 200 


SWP51852LL 


large 

1994-08-20 09:40:09 

900 


LWP30479LL 


large 

1995-04-19 02:31:12 

600 


SWP54463LL 


large 

1995-04-19 01:59:00 

1200 


SWP54465LL 


large 

1995-04-19 08:18:31 

1800 


LWP30490LL 


large 

1995-04-20 01:45:20 

1800 


SWP55690HL 


large 

1995-08-24 18:23:29 

22920 


LWP32551LL 


large 

1996-08-10 20:08:00 

1200 


LWP32552LL 


large 

1996-08-10 21:54:07 

1800 


SWP57774LL 


large 

1996-08-10 19:31:12 

1800 


SWP57775LL 


large 

1996-08-10 20:55:02 

3000 

FOS, HST 

Y3415205T 


073, G190H 

1997-02-03 11:21:13 

180 


Y341520BT 


073, G400H 

1997-02-03 12:25:41 

950 


Y341520CT 


0/3, G570H 

1997-02-03 12:48:27 

1500 


Y341520KT 


0/3, G270H 

1997-02-03 15:38:38 

1650 

STIS, HST 

04NH01030 

52" 

x 0/05, G750M 

1998-03-07 22:58:48 

708 


04NH01040 

52" 

x 0/05, G430L 

1998-03-07 23:15:55 

840 


04NH02030 

52" 

x 0/05, G230LB 

1999-02-05 18:13:25 

1978 


04C562020 

52" 

’ x 0/05, G230L 

1999-04-12 01:22:42 

900 


04C562030 

52" 

x 0/05, G140M 

1999-04-12 01:45:28 

862 


06IH20020 

52" 

x 0/05, G430M 

2001-09-16 02:35:05 

900 


06IH20060 

52" 

x 0/05, G430M 

2001-09-16 03:20:42 

1080 


06IH30020 

52" 

x 0/05, G750M 

2001-09-24 08:55:01 

1080 


06IH30060 

52" 

x 0/05, G750M 

2001-09-24 10:32:46 

1080 

FUSE 

Q3040101000 


LWRS 

2002-04-12 02:00:59 

8240 


Q3040102000 


LWRS 

2002-06-16 22:12:47 

15 268 


U 109320 1000 


LWRS 

2006-06-23 04:02:36 

15 824 


U 1093202000 


LWRS 

2006-06-24 11:12:43 

2862 
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Table 3. Statistics of our model atoms. 


Ion 

NLTE levels 

Lines 

Sample lines 

H 

I 

5 

10 



II 

1 

- 


He 

I 

29 

69 



II 

20 

190 



III 

1 

- 


C 

III 

58 

329 



IV 

54 

295 



V 

1 

0 


N 

III 

34 

129 



IV 

90 

546 



V 

54 

297 



VI 

1 

0 


O 

III 

72 

322 



IV 

38 

173 



V 

53 

255 



VI 

54 

291 



VII 

1 

0 


Ne 

II 

14 

19 



III 

14 

17 



IV 

14 

24 



V 

14 

18 



VI 

1 

0 


Si 

III 

17 

28 



IV 

16 

44 



V 

25 

59 



VI 

1 

0 


P 

IV 

15 

9 



V 

18 

12 



VI 

1 

0 


S 

IV 

20 

25 



V 

23 

33 



VI 

30 

42 



VII 

1 

0 


Fe 

IV 

5 

12 

3 102 371 


V 

7 

25 

3 266247 


VI 

7 

25 

991935 


VII 

1 

0 


Ni 

IV 

5 

14 

2512561 


V 

7 

25 

2766664 


VI 

7 

22 

7408 657 


VII 

1 

0 



Total 

467 

1022 

20048 435 


Notes. For iron and nickel, we used a statistical superline/superlevel 
approach (Rauch & Deetjen 2003) in the subsequent NLTE level popu- 
lation calculations. The original number of the so-called “sample lines” 
from Kurucz’s line list, which are combined to super-lines, is also given. 



Fig. 6. Comparison of the temperature structures of a HHeCNO (blue, 
thin) and the final HHeCNONeSiPSFeNi model (red, thick) with 
T eS = 60 kK and logg = 6.0. 



ax/ A 

Fig. 7. Determination of the Si abundance. Synthetic line profiles cal- 
culated from models with O.lx solar (blue, dashed), solar (green, light 
gray, including ISM absorption lines), and lOx solar (red, dark) Si abun- 
dance are compared with the observed Si IV (left) and Si III (right) line 
profiles. 



Fig. 8. Determination of the P abundance. Synthetic line profiles calcu- 
lated from models with 0.0 lx solar (blue, dashed), O.lx solar (green, 
light gray, including ISM absoiption lines), and solar (red, dark) P abun- 
dance are compared with the observed P V line profiles. 
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Fig. 9. Determination of the Fe abundance and the upper limit for the Ni abundance. Synthetic spectra calculated from models with O.lx solar 
(blue, dashed), solar (green, light) and lOx solar (red, dark) Fe and Ni abundances are compared with the observation. 
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